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We have measured a paramagnetic Meissner effect in Nb — AI2O3 — Nb Josephson junction arrays 
using a scanning SQUfD microscope. The arrays exhibit diamagnetism for some coohng fields 
and paramagnetism for other cooling fields. The measured mean magnetization is always small, 
{$tot — $ext) < O.S^o (in terms of flux per unit cell of the array, where "l>o is the flux quantum) for the 
range of cooling fields investigated ( — 12$o to 12cl?o). We demonstrate that a new model of magnetic 
screening, valid for multiply-connected superconductors, reproduces all of the essential features of 
paramagnetism that we observe and that no exotic mechanism, such as d-wave superconductivity, 
is needed for paramagnetism. 



PACS numbers; 74.25.Ha, 74.2Q.-z, 74.50. -l-r, 75.20.-g 



l^The observation of PME 
i2l known to be s-wave sup- 



A paramagnetic response was first observed in 
BSCCOEl and, subsequently reported in many other high 
Tc materials. 0^13 This response is in striking contrast to 
the diamagnetic response which is a hallmark of super- 
conductivity. It has been argued that the paramagnetic 
Meissner effect (PME) results,from a d-wave symmetry in 
the BSSCQ. m'der parameter 
in niobiumLrLi and aluminum, 
erconductors, questioned the vahdity of the d-wave argu- 
ment. 

There are several proposed theories for the PME, al- 
though no experiment has been yet able to provide con- 
clusive evidence for one or another; and there is disa|gree- 
ment as to the validity of the d-wave explanation.^ In 
high Tc materials one possible explanation is the exis- 
tence of TT-junctions which result from the rpisalignment 
of grains in materials with d-wave symmetry.Q In the con- 
ventional materials, which are known to be s-wave, PME 
has been attributed to different mechanisms causing non- 
equilibrium flux configurations. A nonT^jquilibrium flux 
state may result faps fiux-compression,t3 the similar gi- 
ant vortex stateJi3o pinning raadomness, or a surface 
barrier at the edge of the sample.ll3 

To create a controlled experiment we have chosen 
to look at Josephson-junction arrays using a scanning 
SQUID microscope (SSM), using a technique similar to 
that of Kirtley et a/.li3 The SSM provides spatially re- 
solved magnetization images of the sample. Arrays elim- 
inate many of the complications in studying the underly- 
ing cause of the PME. Since our arrays are lithographi- 
cally fabricated from Nb and AI2O3, we know that there 
are neither 7r-junctions nor significant randomness. Ad- 
ditionally, our sample is uniformly attached to the ther- 
mal bath, eliminating thermal gradients across the sam- 
ple during cooling. We find that arrays can be either 
paramagnetic or diamagnetic depending upon the cool- 



ing field. Additionally, we propose a new model, valid 
for multiply-connected supercondutors, which predicts 
a paramagnetic response and agrees with our scanning 
SQUID measurements. 

The arrays consist of a square array of niobium crosses 
in two layers, with a unit cell size of a = 46/im. Junc- 
tions are formed at the cross overlap. The calculated self- 
inductance of each loop of four junctions is i = 64pII and 
the measured critical current density is Jc = 600A/cm^ 
at 4.2K, with a junction area of 5 x 5/zm^. The entire 
array consists of 30 x 100 junctions. The unit-cell dimen- 
sions comnare to the typical grain size seen in BSCCO 
sampleaJtj which exhibit PME, to which similar fields 
are applied. The array is placed in an SSM sample 
stage,EZl sitting at the end of a sapphire rod of 10mm 
diameter, around which a solenoid has been wound. Be- 
fore any data is collected we first remove any flux trapped 
within the pick upJoop of the SQUID to avoid interaction 
with the samplcEj Additionally, we measure and correct 
for any background field present. 

We cooled our arrays to 4.2K in fixed external mag- 
netic field, and then used the SSM to measure the mag- 
netic flux threading the SQUID from the sample, ^tot, 
as a function of position. Our spatial resolution was lim- 
ited by our SQUID-sample separation distance, which we 
estimate to be between 40 and 60/im. Data were taken 
every 5/im in the x direction and every 50/im in the y 
direction, and converted from flux through the SQUID 
to normalized flux per unit cell. We measure the exter- 
nal flux, $oxt, directly by warming the sample above 
and imaging. From the data we calculate both a local 
normalized magnetization ($tot — 'I'oxt)/*I*o an aver- 
age overall magnetization ($tot — ^'ext}/$o-lij We-svere 
unable to resolve individual unit cells in the array,c2l but 
we were still able to resolve large scale flux distributions. 

Our measurement and analysis proceeded in a fash- 



ion similar to Ref. We first make an image of the 
magnetization of the array (a typical image is shown in 
Fig. |l|(a)). We then make a magnetization histogram 
and compute the mean magnetization. Fig. ^b) shows 
the magnetization histogram generated from the image 
in Fig. H(a). The histogram is computed by looking at 
the lower third of the array and cptinting the number 
of points at a given magnetization.Eil We find that the 
array may be either diamagnetic or paramagnetic, but 
surprisingly, appears to be preferentially paramagnetic. 
Additionally, we always observe diamagnetic screening 
currents around the outside edges of the sample. The ar- 
ray shown in Fig. |l|(a) was cooled in an external field of 
4.8$o and clearly is paramagnetic {i.e. it is darker than 
the background). We find a mean paramagnetic response 
of ($tot - $ext) = (0.063 ±0.005)$o for this cooling field. 
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FIG. 1. (a) Representative paramagnetic magnetization 
image for the array, with an external cooling field of 4.8"I>o- 
The color scale runs from — 0.6$o (white) to 0.6$o (black), 
(b) Magnetization histogram generated from the paramag- 
netic image. 

The data plotted in Fig. |^ show the mean magnetiza- 
tion versus external cooling field. Each data point repre- 
sents a single field-cooled average measurement. The er- 
ror in the mean values plotted, such as discussed above, 
is less than the size of each data point as drawn. The 
sign of the magnetization is reproducible, and the mag- 
nitude is reproducible to within a factor of two when the 
sample is warmed above lOK and recooled in the same 
field, indicating that our external field is consistent to 
within iO.OS^o The external field has a spatial variation 
of ±0.5$Oj which is limited by the gradient of the back- 
ground field in our probe. Although there is some varia- 
tion from one cooling to the next in the actual magnetic 
images for the same external field, the paramagnetism is 
remarkably reproducible. 
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FIG. 2. 



-0.C5- 
-0.10- 

-0.15- 

♦ 

-0 20- 

Measured mean magnetization of the array 



($tot - $cxt)/"I'o, plotted versus cooling field, $oxt/"I>o- The 
error for each data point is less than the size of the plotted 
symbols. 

We may consider a single loop of the array as a jumping 
off point for modeling the entire array's magnetization. 
Originally used to describe the observed magnetiaation of 
a single loop with a single Josephson-junctionJHj Sigrist 
and Ricecl later used this model to argue that 7r-junctions 
could cause PME in d-wave superconductors. The single- 
loop model gives a relationship between the external flux 
applied to a four-junction loop and the total measured 
flux 
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(1) 



Here Ic is the critical current of a single junction. For 
all values of > -, <i>tot is a multi- valued function 
of <i>oxt depending on an integer n. We note here that 
Ref. ^ used a model with a slightly greater than one 
to demonstrate the absence of potential paramagnetism 
for a single loop with a single junction. Had they used a 
value of ^ 1 they would have noticed the potential 
for paramagnetism even in the absence of external fiux, 
although such a situation is not energetically favorable. 
For a four junction loop, Pi = 27r-^ = 30 allows for 
the possibility of paramagnetism (or diamagnetism) at 
$cxt = 0. 

Computing the Gibbs free energy for the various possi- 
ble solutions to (|^) we find that for m < < the 

single loop will be diamagnetic and for m -f- i < < 
m + 1 will be paramagnetic, where m is any integer. If 
one applies a positive offset fiux <i>off to the four-junction 
loop {i.e. the total applied flux is <I>cxt + 'i'off) the so- 
lutions $tot — ^cxt vs. $oxt become more paramagnetic 
and conversely if $off < the solutions become more dia- 
magnetic. If, e. 5., a simple diamagnetic screening current 
flows around the outside edge of the superconductor (neg- 
ative offset), the measured magnetization will be more 
diamagnetic. 

A conventional superconductor is a perfect diamagnet: 
An applied fiux causes screening currents to fiow around 



the outside. In an ideal array of Josephson junctions, 
one may decompose this edge current, /, into loop cur- 
rents of magnitude /, flowing around each plaquette, as 
shown in Fig. ||(a) . In a real array, this is not likely to be 
the complete picture. The currents will flow only within 
one penetration depth in the superconductor. For our 
array, the penetration depth is 900A while the wire size 
is lO^m so that the interior loop currents do not cancel. 
This means that the array must create an energy ^LI^ 
for each plaquette in order to generate a screening cur- 
rent / for the array. In an overlap geometry, the currents 
through the Josephson-junctions are uniform, so the cur- 
rents through the junctions still cancel and we do not 
have any increase in the total Josephson energy, Ej. For 
an array oi N xN junctions Fig. ^a) costs a total energy 
of 
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(4) 



Ea = -LI^N^ + 4iV(l - cos-/)Ej 
where 7 is the gauge invariant phase difference. 
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ing current around a multiply connected superconducting 
sample: (a) The conventional picture, where each loop of the 
sample generates a current, (b) Our new picture, where only 
the exterior plaquettes create a screening current. This cre- 
ates a diamagnetic current on the outside of the sample, and 
a paramagnetic current just inside the sample. 

We propose that the array actually screens in a much 
different fashion. The array may generate a diamagnetic 
screening current around the outside perimeter if each ex- 
terior plaquette maintains a screening current, as shown 
schematically in Fig. ||(b). In addition to the diamagnetic 
screening current, we then also have a paramagnetic cur- 
rent, of the same magnitude, flowing just inside the array. 
The array must generate the ^LI^ energy only for each 
exterior plaquette, so the energy cost to the array in this 



Et = -LI^{AN) + {4A^ + 4(iV 



l)}(l-cos7)Sj. (3) 



It is clear from (H) and (|) that Ea oc N'^ while Eb (x N 
so that for a large enough array, the conventional picture 
of screening will not be energetically favorable. Fig. ||(b) 
has lower energy than Fig. }Q{a) when, assuming 7 is pro- 
portional to $tot, 



For our array the criteria is easily satisfied and we expect 
screening to take place as in Fig. ||(b). 

In Fig. ^(a) we can treat the outer ring of Josephson- 
junctions as a single loop of 4iV junctions. This loop will 
screen up to — 4.8 flux quanta from the interior of the 
array, when Ic is flowing around the exterior of the array. 
In Fig. H(b) each exterior plaquette can screen up to 
flux quanta from itself. This is a much larger absolu1;e 
fleld than in the conventional case and the difference can 
be easily observed experimentally. We have observed, 
that flux will not penetrate the array until the external 
flux exceeds at least 4.0$o per unit cell. Additionally, 
ac susceptibility work demonstrated that little effect was 
seen until the amplitude of the|-gLriving ac flux exceeded 
4.6$o per unit cell of the arrayH 

The interior plaquettes see the external field, a diamag- 
netic flux from the extreme outside loop of the array, and 
a paramagnetic fiux from just inside the array. Because 
the loop generating the paramagnetic flux is closer, the 
interior loops will see a net paramagnetic flux in addi- 
tion to the external field. Considered individually, the 
interior plaquettes obey a modified Eq. (1h), 
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(5) 



where ^scr is the flux seen by the inner plaquette due 
to the screening currents produced by the exterior pla- 
quettes. When $scr > the interior plaquette sees a 
paramagnetic offset and its magnetization will be more 
paramagnetic. This causes the array to be paramagnetic 
more often than diamagnetic in agreement with experi- 
ment. 

If we consider only the mutual inductance between a 
test plaquette in the interior of the array and the exterior 
plaquettes which provide the screening, we can make a 
crude estimate of the paramagnetic response. The maxi- 
mum current that the junctions can sustain is the critical 
current, 1^, this sets the upper bound for the magneti- 
zation offset. We have numerically computed the flux 
induced from the screening currents in Fig. ^(b) onto 
the interior plaquettes for our 30 x 100 junction array. 
When the screening currents are Ic, the maximum flux 
is induced, and the minimum <I>scr we compute is O.lS^o 
per unit cell, near the center of the array. This means 
that the center plaquette should have its magnetization 
shifted upward by O.IS^q. Plaquettes nearer to the exte- 
rior see $scr values as large as several $o- Additionally, 
we compute an average flux induced into the interior of 
the array to be <I>scr = 0.277<i>o per unit cell. Considering 
the crudeness of the model, this is remarkable agreement 
with our measured magnetization (see Fig. ^). 

We expect that the screening current will increase 
monotonically with the external flux, up to Ic, and there- 
fore expect that the paramagnetic offset in the interior 
plaquette magnetization will increase as the external flux 



increases. This is also consistent with our observations 
in Fig. |. 

The model described here does not account for many 
features of the real array. We have not accounted for 
variations in the Josephson-junction parameters and we 
have ignored the mutu^induction between interior pla- 
quettes. Phillips et alM showed the importance of in- 
cluding the mutual inductance matrix when computing 
array magnetization. Chandran has, in the overdamped 
junction limit, used a model including the full induc- 
tance matrix and computed— a paramagnetic moment in 
a Josephson-junction array.c3 Unfortunately, our arrays 
are highly undcrdampcd, which makes a model utilizing 
the fuU-inductarice and dynamics much more computa- 
tionally difficult .Ej The simplicity of our model does not 
describe all of the features that we observe in our array. 
The observed flux pattern is very complicated, with both 
paramagnetic and diamagnetic regions. 

We have shown that in a multiply connected supercon- 
ductor such as our array, that the lowest energy configu- 
ration that provides diamagnetic screening will also tend 
to make the interior of the array paramagnetic. Further- 
more, the predictions of this model are consistent with 
our measurements, that the array is more often para- 
magnetic than diamagnetic. We believe that the param- 
agnetism in Josephson-junction arrays is best described 
by the effects detailed above. 
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